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CN : Abstract 

f^ ' The modification of mesonic observables in a hot medium is analyzed as a tool to investigate the 

^D ■ restoration of chiral and axial symmetries in the context of the Polyakov-loop extended Nambu- 

- ^^ , Jona-Lasinio model. The results of the extended model lead to the conclusion that the effects 

X 

$H ' of the Polyakov loop are fundamental for reproducing lattice findings. In particular, the partial 

restoration of the chiral symmetry is faster in the PNJL model than in the NJL one, and it is 
responsible for several effects: the meson-quark coupling constants show a remarkable difference 
in both models, there is a faster tendency to recover the Okubo-Zweig-Iizuka rule, and finally the 
topological susceptibility nicely reproduces the lattice results around T/Tc «:: 1.0. 
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I. INTRODUCTION 

In recent years several studies have been developed, which are concerned with the prop- 
erties of matter under extreme conditions of density and/or temperature: the restoration of 
symmetries (e.g., the chiral symmetry) and the phenomenon of deconfinement, which might 
be achieved in ultrarelativistic heavy-ion collisions or in the interior of neutron stars, de- 
serves special attention. Properties of hadrons, in particular, of mesons, propagating in a hot 



m. 



or dense medium can shed light on the occurrence of the expected phenomena [l|, |2|, 

For example, a criterion to identify an effective restoration of chiral (axial) symmetry is to 

look for the degeneracy of the respective chiral partners |6| . 

The study of meson properties, in the SU/(2) sector, around the critical region where the 
phase transition takes place, was performed in Ref. J7|, in the framework of the modified 



Nambu-Jona- 
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jasinio mode 
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including the Polyakov loop (the so-called PNJL model) 
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161 . llTl . Il8l | . In the PNJL model, quarks are simultaneously coupled to 



the chiral condensate and to the Polyakov loop, so the model incorporates features of both 
chiral and Z^r^ symmetry breaking [71]. The coupling to the Polyakov loop is fundamental 
for reproducing lattice results on QCD thermodynamical quantities pLSJ, since it produces 
a suppression of the unphysical colored states (one or two quark states), which should not 
contribute to the thermodynamics below the critical temperature. 

In this paper, we intend to extend the investigation of light scalar and pseudoscalar 
mesons at finite temperature, generalizing previous works to the SU/(3) sector. In particular, 
it will be interesting to compare the properties (e.g., the masses) of the scalar mesons 
(o", /o, ao, and Kq) with those of the pseudoscalar nonet {ri,ri',Tr^, and K), which can be 
considered as chiral partners of the former. We focus our attention on the role of the 
Polyakov loop in determining the behavior with the temperature of these mesons. 

An interesting open question we wish to address is whether both chiral SU(A^/)®SU(A^/) 
and axial 11^(1) symmetries are restored and which observables could carry information 
about these restorations. Moreover, we shall investigate the role of the Ua(1) anomaly 
which, as is well known, is responsible for the flavor mixing effect that removes the degen- 
eracy among several mesons. The 11^(1) symmetry does not exist at the quantum level 
being explicitly broken by the axial anomaly ^] which, in turn, can be described at the 



semiclassical level by instantons 20|. 



The flavor mixing induced by the presence of the axial anomaly causes a violation of the 
Okubo-Zweig-Iizuka rule, both for scalar and pseudoscalar mesons, hence the restoration of 
axial symmetry should have relevant consequences on the phenomenology of the mesonic 
mixing angles as well as on the topological susceptibility. Also in this context the addition 
of the Polyakov loop might influence the tendency toward the recovery of the pseudoscalar 
and scalar mixing angles, already evaluated within the pure NJL model. 

The paper is organized as follows: in Sec. II, we present the model Lagrangian of PNJL 
in SU/(3) including the 't Hooft interaction term; specific subsections are dedicated to the 
gap equations, and to the properties (masses and mixing angles) of the pseudoscalar and 
scalar meson nonets. In Sec. Ill, we show our results, starting with a discussion on the 
characteristic temperatures and the role played by the strange components; then we display 
the mesonic masses, the meson-quark coupling constants, and the topological susceptibility. 
Our conclusions are reported in Sec. IV. 

II. MODEL AND FORMALISM 

A. The PNJL model with three flavors 

We perform our calculations in the framework of an extended SU/(3) PNJL Lagrangian, 
which includes the 't Hooft instanton induced interaction term that breaks the Uyi(l) sym- 
metry; moreover quarks are coupled to a (spatially constant) temporal background gauge 



field representing the Polyakov loop 21 
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+ go {det [q (1 + 75) q] + det [q (1 - 75) q]} 

-UmA],^A];T). (1) 

Here, q = {u, d, s) is the quark field with three flavors {Nf = 3) and three colors {Nc = 3), 
m = dmg{mu,md,ms) is the current quark mass matrix, and A'' are the flavor SUj(3) Gell- 
Mann matrices (a = 0,1,..., 8), with A° = \/|l- The covariant derivative is defined as 
£)M — Qfj. _ j^j^i-i^ with A^^ = 5qA^ (Polyakov gauge); in Euclidean notation y4° = —iA^. The 
strong coupling constant G strong is absorbed in the definition of A^(x) = G strong-^ai^) ^ ^ 
where A^ is the (SUc(3)) gauge field and A^ are the (color) Gell-Mann matrices. 



The Polyakov loop field $ appearing in the potential term of (ITD is related to the gauge 



field through the gauge covariant average of the Polyakov line 10 



13| 



$(x) = ((/(f))) = — Tr,((L(f))) 



Nr_ 



where 



L(x) = V exp 



I I dTA4{x,T) 




(2) 



(3) 



The Polyakov loop is an order parameter for the restoration of the Z3 (the center of SUc(3)) 
symmetry of QCD and is related to the deconfinement phase transition: Z3 is broken in the 



25 



26| 



deconfined phase ($ ^ 1) and restored in the confined one ($ ^ 0) [2J, 

Here, it is important to make some remarks about the applicability of the PNJL model. 
Beyond the chiral pointlike coupling between quarks, in the PNJL model, the gluon dynamics 
is reduced to a simple static background field representing the Polyakov loop (see details 
in Refs. [71, ll3|]). This scenario cannot be expected to work outside a limited range of 
temperatures. Indeed, at large temperatures it is expected that transverse gluons start to 
be thermodynamically active degrees of freedorn, but they are not taken into account in 
the PNJL model. Since, as concluded in Ref. 27|, transverse gluons start to contribute 
significantly for T > 2.5 T^, where T^ is the deconfinement temperature, we can assume that 
the range of applicability of the model is roughly limited to T < (2 — 3)Tc. 

Concerning the effective potential for the (complex) $ field, there exist in the literature 



different choices 
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23|: we decided to adopt the one proposed in Ref. |15|] [see Eq. H], 
which is known to give sensible results |l5l . [l6|. In particular, it reproduces, at the mean 
field level, results obtained in lattice calculations. The potential reads 



U (^ ^-T) n (T) - 



where 



y4 



a (T) = ao + ai ( Y ) + «2 ( ^ ) and b{T) 



b.i^ 



(4) 



(5) 



We notice that in the mean field approximation the Polyakov loop field ^{x) simply 
coincides with its expectation value $ =const., which minimizes the potential (jl]). The 
parameters of the effective potential U are given in Table [B These parameters have been 
fixed in order to reproduce the lattice data for the expectation value of the Polyakov loop 
and QCD thermodynamics in the pure gauge sector 28|, |29 |. 



ao 


ai 


0-2 


^3 


3.51 


-2.47 


15.2 


-1.75 



TABLE I: Parameters for the effective potential in the pure gauge sector (Eq. ([1])). 



The parameter Tq is the critical temperature for the deconfinement phase transition within 
a pure gauge approach: it was fixed to 270 MeV, according to lattice findings. Different 
criteria for fixing Tq are available in the literature, like in Ref. [181], where an explicit Nf 
dependence of Tq is presented by using renormalization group arguments. Besides, we notice 
that the Polyakov loop computed on the lattice with (2-1-1) flavors and with fairly realistic 
quark masses is very similar to the SU/(2) case 29|. Hence, we choose to keep for the 
effective potential U ($, $; T) the same parameters which were used in SUj(2) PNJL 15| . 
including Tq. 

Since one of the purposes of the present paper is to estimate the effect of the coupling 
of the Polyakov loop to quarks with NJL 4-point interactions, we choose to compare the 
NJL and PNJL models on a relative temperature scale T/T^, where T^ is a characteris- 
tic temperature that can be derived in both models (here, the chiral crossover location). 
This choice is justified along the lines of the Ginsburg-Landau effective theory, where char- 
acteristic temperatures cannot be absolutely compared between two models of the same 
universality class. Besides, as noticed in Refs. [Tl, [l8|] the Tq dependence of the results is 
mild. In the present context the physical outcomes are not dramatically modified when one 
changes Tq^. Hence the choice Tq = 270 MeV appears to be preferable here, since it ensures 
an almost exact coincidence between chiral crossover and deconfinement at zero chemical 
potential, as observed in lattice calculations. 

With the present choice of the parameters, $ and $ are always lower than one in the 
pure gauge sector. In any case, in the range of applicability of our model (T < 2.5 Tc), there 
is a good agreement between our results and the lattice data for $. 

Let us anticipate that, at T = 0, it can be shown that the minimization of the grand 
potential leads to $ = $ = 0. So, the quark sector decouples from the gauge one, and the 
model is fixed by the coupling constants gs, go, the cutoff parameter A, which regularizes 



^ The choice of Tq certainly deserves investigations beyond the scope of this paper. 



the divergent integrals, and the current quark masses rrii. The parameter set used here is 
mu = md = 5.5 MeV, m, = 140.7 MeV, gsA^ = 3.67, goA^ - 12.36 and A = 602.3 MeV (for 
details see Ref. |30|). 



B. Gap equations 

From the Lagrangian ([T]) in the mean field approximation it is straightforward (see 
Ref. [311) to obtain effective quark masses (the gap equations) given by 



Mi = mi- 2gg {qiQi) - 1g^ {q^q^) {q^qt) 



(6) 



where the quark condensates {qiqi), with i,i,k = u,d,s (to be fixed in cyclic order), have 
to be determined in a self-consistent way. The last term on the right-hand side derives 
from the determinantal piece in the Lagrangian, which clearly induces a fiavor mixing in 
the constituent quark masses Mj. So, taking already into account Eq. ([6]), the PNJL grand 
canonical potential density in the SU/(3) sector can be written as 

Q = f]($, $, M,; T,fi)=U ($, $, T)+g^ ^ (g,g,)' + Ag^ {q^q^) {q,q,) (g, g,) 

{i=u,d,s} 

-2^^ E /7^^'-2T Yl I PrA4{Ed + zliE^)) , (7) 



{i=u,d,s} ' 



(27r) 



{i=u,d,s} ' 

where we have defined E^ = Ei ^ ft, the upper sign applying for fermions and the lower 
sign for antifermions; Ei is the quasiparticle energy for the quark i: Ei = a/p^ + M^; finally, 
2;J and z^ are the partition function densities (with the usual notation f3 = 1/T) 



z+{Ei) = Tr.ln 
z^{Ei) = Trcln 



1 + Lte-^^^ 



(+) 



-m 



(+) 



A-y 



In <i 1 + 3 ( $ + <l>e '^"> ) e 



-m 



(+) 



+ e 



-313 El 



(+) 



(8) 
1 + Le-f'^r^ In |l + 3 ($ + l>e-^^r') e"^^'"' + e'^^^'"' } . (9) 

It was shown in Ref. [Tj that all calculations in the NJL model can be generalized to the 
PNJL one by introducing the modified Fermi-Dirac distribution functions for particles and 
antiparticles: 



/$ {Ei 



exp{z^{Ei)} 

$e-/^^l'' + 2$e-^^^'"' + e-=^^^'"' 
exp{z^{Ei)} 



(10) 

(11) 



To obtain the mean field equations we must search for the minima of the thermodynamical 
potential density ([7]) with respect to {qiqi) {i = u,d,s), ^, and $. In fact, by minimizing Q 
with respect to {(jiqi), we obtain the equations for the quark condensates 

fe) = -2nJ ^^[l - ft\E.^) - fi-\E,)] . (12) 

J A [ZTr) J^i 

Furthermore, minimization of fl with respect to $ and $ provides, respectively, the two 
additional mean field equations [7| 



a{T) - b{T) [$ - 2<l>2 + $2$] 




1 - 6$$ + 4(<l>3 + $3) - 3(<l'<l')2 



+ ^^77TTt1 (13) 



$-6 



(27r)^ yexp{4(Ei)} exp{z^{Ei)} 
h{T) [$ - 2l>2 + $$2] 



1 - 6$$ + 4(<l>3 + $3) - 3($$)2 
e '^ » e 



"^^ {.=M,.} -^^ (2^ Vexp{4(E,)} + exp{^^(E,)} ) ^^^^ 



This general formalism, presented here for completeness in the grand canonical approach, is 
going to be employed in the present work with n = 0. Under this condition, the constraint 
$ = $ holds. 

C. Pseudoscalar and scalar meson nonets 

To calculate the meson mass spectrum, we use the same procedure described in detail in 
Ref. |6|. The model Lagrangian ([1]) can be put in a form suitable for the usual bosonization 
procedure, after reducing the six-quark interaction term in order to obtain an effective 
four-quark interaction. This can be achieved, e.g., by shifting the operator {qX"'q) — > 



{qX"'q) + (^A"g), where (...) is the vacuum expectation value [3l|| and by Wick contracting 
one of the bilinears {qX"'q) in the resulting six-quarks interaction term. The following effective 
quark Lagrangian is thus obtained: 

Ceff = q{i-f^d^ - rh)q 

+ ^5a4(gA"g)(gA^)]+ ^P„5[(g^75A'^g)(g^75A^)], (15) 



where the projectors Sab , Pat are of the form 

Sab = gsSab + guDabc {q^q) , (16) 

Pab = gsSab - goDabc {q>^''q) ■ (17) 

The constants Dabc coincide with the SU/(3) structure constants for a,b,c = (1,2, ... ,8), 
while Doab = —y^Sab and Dqoo = y §• The bosonization procedure is then reahzed by 
integrating out the quark fields in the functional generator associated with the Lagrangian 
flTSjl complemented by the coupling to scalar and pseudoscalar bosonic fields. With this 
procedure the natural degrees of freedom of low-energy QCD in the mesonic sector emerge 
in the resulting effective bosonic action. By expanding the latter up to second order in 
the meson fields, one obtains the meson propagators, from which the masses, meson-quark 
coupling constants and meson decay constants /m can be evaluated according to the usual 



definitions (see, for example, Eqs. (21) and (22) of Ref. 32l|). 



For example, we obtain the inverse propagator for the pion as 

D-\P) = 1-P^U:^{P), (18) 

where 

P7r = gs + go {qsqs) , (19) 

and n^j(P) is the polarization operator, which in momentum space has the form 

Kb{P) = ^N, I ^,trn [S.ip)iX''Utl,)S,{p + P)(A^),.(^75)] , (20) 

with tr^) being now the trace over Dirac matrices. For details concerning the calculations for 







the other meson propagators, see Ref. [6[. At variance with Ref. [6|, here the introduction 
of the Polyakov loop implies, obviously, the use of the modified Fermi functions /$ in the 
calculation. 

The inclusion of the 't Hooft interaction in the PNJL/NJL model allows for flavor mixing; 
however, by imposing SU/(2) flavor symmetry (namely, rriu = rrid) the off-diagonal coupling 
strengths that mix the 7r° with rj and 1]' vanish; hence, charged and neutral pions become 
degenerate in mass, as well as the neutral and charged kaons. 

It should be noticed that flavor mixing somewhat entangles the calculation of the r/ and 



rj' masses and couplings, since it gives rise to a P^-dependent mixing angle 9p{P'^) p, l30|, 



which in our scheme is defined as follows: 

= O(0p)f*Uf"^*''-""*''U* |. (21) 

7]' j \ ^0 / \ ^^^^^ ^°^^^ / \ '^^ 

In the above, r] and 77' stand for the corresponding physical fields, while 773 and 770 are the 
mathematical objects transforming as octet and singlet states of the SU/(3) pseudoscalar 
meson nonet, respectively. By using a standardprocedure we get the equation for the mixing 
angles and for the inverse meson propagators 6|. 

As shown elsewhere, in the framework of the NJL model [g, |30(], the mixing angle 9p^ 
between the components 770 and 7^8 is P^-dependent; hence, one gets different mixing angles 
at P^ = M^ or P^ = M3. In the present paper we only discuss the mixing angle at 
P^ = M^; nevertheless, we checked that the behavior of the mixing angle at P^ = M3 gives 
qualitatively similar information. 

The same technique used for the pseudoscalar sector can now be directly applied to the 
scalar resonances. We deal here with nine scalar resonances: three ao's, which are the scalar 
partners of the pions, four k's, being the scalar partners of the kaons, and the a and /o, 
which are similarly associated with the rj and 7/'. As in the pseudoscalar case, we impose the 
SU/(2) flavor symmetry, so we have no mixing between the a and /o and the neutral Oq. 

As is well known, when the mass of the meson exceeds the sum of the masses of its 
constituent quarks, the meson can decay in its quark-ant iquark pair, hence becoming a 
resonant state. So, in order to calculate the mass of the resonance, the imaginary part of 
the propagator must be taken into account as well and, following the standard approximation 



described in Ref. 33|], one can obtain the mass and the corresponding decay width. 



III. NUMERICAL RESULTS 



A. NJL vs PNJL: Characteristic temperatures 

We start ou. analysis by identifying the chatactetistic^ten^petatutes that separate the 
different thermodynamic phases in PNJL and NJL models p]. In the framework of the PNJL 
model the critical temperature related to the "deconfinement"^ phase transition is T* and 



^ The terminology "deconfinenient" in our model is used to designate the transition between $ ~ and 
$ ~ 1 (see Ref. (7| for a discussion of this point). 
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FIG. 1: Left panel: comparison of the quark masses in the PNJL (sohd hnes) and NJL (dashed 
hnes) models as functions of the temperature; the Polyakov loop crossover is also shown. Right 
panel: the same as before as a function of the reduced temperature T/Tc- 



corresponds to the $ crossover location (FigH]). The chiral phase transition characteristic 
temperatures, T^, are given in both models by the inflexion points (chiral crossover) of the 
chiral condensates {qiQi). Since in both situations chiral symmetry is explicitly broken by the 
presence of nonzero current quark mass terms, chiral symmetry is realized through parity 
doubling rather than by massless quarks. The effective chiral symmetry restoration [in the 
SU/(2) sector] is signaled by the degeneracy of the chiral partners (vr, a) and {rj, oq) or, 
strictly speaking, by the merging of their spectral functions p, LZj. 

In order to compare the NJL and PNJL results, it is useful to follow the evolution of 
the observables as functions of the temperature, expressed both in physical units (MeV) 
as well as rescaled in units of a characteristic temperature. For the latter we choose the 
corresponding T^ in NJL and PNJL, and most figures will be shown versus T/T^, which 
allows a better understanding of the relevant differences between NJL and PNJL models. 
Indeed, we are not interested in discussing absolute quantities but rather in comparing the 
behavior, below and above T^, of the mesonic properties; the key point under investigation 
is the effective restoration of symmetries, as well as the influence of the Polyakov loop on 
the phase transition. 

In Table [TTl we quote the characteristic temperatures: the effective chiral symmetry 
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Model 


T<? [MeV] 


T* [MeV] 


TLott [MeV] 


niott [MeV] 


PNJL 


250 


233 


267 


237 


NJL 


196 


— 


212 


160 



TABLE II: Characteristic temperatures in the NJL and PNJL models at Tq = 270 MeV and zero 
chemical potential. 

restoration temperature, the deconfinement temperature, and the Mott temperatures for 
the pion and the sigma both in PNJL and NJL models. We remind the reader that the 
Mott transition is associated with the composite nature of the mesons: at the Mott temper- 
ature the decay of a meson into a qq pair becomes energetically favorable. 

As already noticed in the SUj(2) PNJL model 13|, T* differs by only a few MeV's from 
T^. In SU/(3) this difference is 17 MeV. From Table [TTl we also see that the character- 
istic temperatures obtained in PNJL are much larger than the lattice result for the chi- 
ral/deconfinement phase transition in (2+1) flavors QCD [T^ ~ 170 MeV). It was pointed 
out in Ref. 7j] that the difference between T* and T^ is due to the choice of the regular- 
ization procedure; we apply here the three-dimensional momentum cutoff to both the zero 
and the finite temperature contributions. Notice also that a different type of regularization 
can lower T^ 3J]. Another important aspect contributing to the present high value of T^ is 
the fact that we do not rescale the parameter Tq to a smaller value after introducing quarks 
in th e sy stem. Indeed, we checked that using the lower value Tq = 187 MeV suggested in 



Ref. p^, smaller characteristic temperatures, T* and T^ are obtained. However, we prefer 
to adopt the higher Tq value since it gives a smaller difference between the critical points of 
the two transitions. Nevertheless, once we are interested in the general properties of mesons, 
the absolute value of the critical temperature is not so relevant: indeed, these properties 
are independent of the specific value of T^ and a different value of Tq does not change the 
conclusions. 



B. The strangeness in PNJL 



as 



In Fig. [T] we plot the masses of the strange and nonstrange quarks and the Polyakov loop 
functions of the temperature. At temperatures around T^, in both models, the mass of 



11 



the light quarks drops to the current quark mass, indicating a smooth crossover from the 
chirally broken to an approximate chirally symmetric phase: we have a partial restoration 
of chiral symmetry. This dropping is more pronounced in the PNJL model than in the NJL 
one, indicating that the transition toward the partial chirally restored phase is faster in the 
former. The strange quark mass shows a behavior very similar to the one of the nonstrange 
quarks, with a significant decrease above T^\ however, its mass is still far away from the 
strange current quark mass. One can say that chiral symmetry shows a slow tendency to 
get restored in the strange sector even if this tendency is faster in the PNJL model. As 
in the NJL model, for what concerns the strange sector |6| and since rriu = rrid < rris, the 
(sub)group SU(2)(S)SU(2) is a much better symmetry of the Lagrangian ([T]). 

Nevertheless, the fact that the masses of the quarks drop faster around T^ in the PNJL 
model is important for the mesonic properties (for example it could modify the survival of 
mesonic bound states in the plasma phase). Moreover, due to the strange components of 
some mesons, the behavior of the strange quark mass (modified by the Polyakov loop) is 
important for their properties, as well as for other observables related to the axial anomaly 



fas noticed in 



3J] the topological susceptibility is strongly influenced by the strange sector) . 



C. Mesonic masses and mixing angles 

In Fig. [2] we plot the masses of the pseudoscalar mesons and of the respective scalar 
chiral partners as functions of the reduced temperature T/T^. The first evidence emerging 
from these figures is that the behavior of the mesonic masses in PNJL looks, qualitatively. 



very similar to the corresponding one in NJL 30|, l35|]. Hence, the results in the two models 



will deserve, in general, similar conclusions, although there are some quantitative differences 
with a non-trivial meaning that will be commented on below. 

In the upper panel of Fig. [2] we plot the masses of the scalar and pseudoscalar mesons 
(cr, ao, vr, rj) (upper panel), /o , V (middle panel) and finally the masses for the K meson 
and its chiral partner n (lower panel), both in the NJL and PNJL models. The lower limits 
of the continua pertaining to each meson are also shown. In fact, the continuum starts at 
the crossing of the vr, a, and rj lines with the quark threshold 2Mu (upper panel), of r/' with 
2Ms (middle panel) and of the K line with Mu + M^. 

We will now analyze the general behavior of the mesons and mixing angles, starting by 
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FIG. 2: Comparison of the pseudoscalar and scalar mesons masses in the PNJL (thick lines) and 
NJL (thin lines) models as functions of the reduced temperature T/Tc ■ In the upper panel the 
oo (dashed line), a (dotted-dashed line), r] and vr (continuous lines) are shown, together with the 
2Mu mass (dotted lines). In the middle panel /o (dashed line) and ??' (continuous line) are shown 
and compared with the 2Ms mass (dotted line). In the lower panel the n (dashed line) and K 
(continuous line) masses are compared with M„ + Mg (dotted line) . 

emphasizing what both models have in common. Concerning the pseudoscalar mesons, it is 
found that they are bound states at low temperature (with the exception of the rj' meson, 
which is always above the continuum ui^ = 2M„), but at the respective Mott temperatures 
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they become unbound (see Table [TTl); as usual, for the vr and K mesons this occurs at 
approximately the same temperature in both models {Tj^j^^^ = 266 (210) MeV). 

The a is the only scalar meson that can be considered as a true (slightly) bound state 
for small temperatures (the other mesons being always resonant states) and turns into a 
resonance at the corresponding Mott temperature (see Table [IT|) . 

Concerning the pseudoscalar mixing angle 9p, it is found that, as the temperature in- 
creases, it approaches the ideal value 6p = —54.736°, although never reaching it (see Fig. [3]). 
As a consequence, the quark content of the mesons i] and t]' changes remarkably, although 
a small percentage of mixing always remains: the t] eventually becomes almost nonstrange, 
while the opposite happens to t]' p, |35| . The scalar mixing angle 9s exhibits a similar ten- 
dency, the ideal mixing angle ds = 35.264° never being reached in the range of temperatures 
studied here (see Fig. |3]); hence, the strange component of the a meson decreases but never 
vanishes, and /o becomes almost purely strange. 

Our main concern in this paper are the modifications introduced by the Polyakov loor 



on the results obtained in the pure NJL model, which are not new (see for instance p, |30|, 



34 



351]); yet, some comments are in order concerning the behavior of the mixing angles, 
keeping also in mind some recent contributions to this subject in the framework of other 
models 



43 



4J] . The mixing angles are very sensitive to the medium effects, in particular to 
the influence of the medium on the strange quark mass: this might also explain why some 
aspects of the in-medium behavior of 6s and 6p are not the same in different models or 
situations. 

First of all it should be noted that the mixing angles depend on the masses of the mesons 
and, for the sake of illustration, the angles plotted here are 9s, depending on the mass of 
the a meson and 6p, depending on the mass of the rj meson. Since the a meson has a small 
strangeness component, its behavior is essentially driven by the decrease of the nonstrange 
quark mass; on the contrary, the t] has a stronger strangeness component, and its behavior 
is affected not only by the fast decrease of the nonstrange quark mass, but also by the slow 
decrease of the strange quark mass. Consequently, the mass of the a decreases more rapidly 
than the one of the rj and, as a result, 6s gets closer than 6p to the respective ideal value, as 
can be seen in Fig. [31 Notice that there is even a slight increase of 6p at about T ~ 1.75 T^, 
the temperature at which the rj meson enters into the strange quark continuum (m^ > 2Ms). 

The evolution of the strangeness content of r], rj' determines which one will become non- 
14 



strange, hence behaving as the chiral axial partner of the vr. For example, the finite tem- 
perature behavior of Op leads to the identification of the t] as the chiral axial partner of 



the vr in Refs. 



40 



4l| and 



yy 



34 



3ut the opposite is found in Refs. J42|, |43|. An 



interesting situation was reported in Refs. 



351 ]. for neutron matter in (3 equilibrium. 



where the strange quark mass decreases more rapidly than in symmetric nuclear matter. 
As a consequence, the pseudoscalar mixing angle changes sign and approaches the positive 
ideal value, while t] and r]' exchange identities, the rj' becoming nonstrange and exhibiting 
a tendency to degenerate with the pion. A level crossing of the pseudoscalar mixing angle 
and exchange of identities of 77, r/ was also found in Ref . 4l(| , and a similar effect for the 
scalar angle is reported in Ref. 431]. 

Recently, a study of the mixing angles in the framework of a Schwinger- Dyson ap- 
proach 4^, using a separable interaction, established a connection between the behavior 
of Op and the fastness of restoration of the axial symmetry, characterized by the drop of 
the topological susceptibility. Since in the NJL model the topological susceptibility is pro- 
portional to the strange quark condensate, this result is compatible with our remark, which 
relates the behavior of the mixing angles to the strange quark mass evolution. We observe 
that in the PNJL model there is a faster restoration of chiral symmetry, both in the non- 
strange and strange sectors, leading to a modification of the in-medium behavior of the 
mixing angles, meson masses, coupling constants, and topological susceptibility, as we will 
discuss below. However, the modification induced by the Polyakov loop is limited to a range 
of temperatures around the critical one (~ 0.75 T^'^ — 1.5 T^^), and it is not strong enough to 
alter the sign of the mixing angles. Having in mind the lack of experimental information on 
the behavior of the observables discussed here, a comparative study of 0$ and Op in different 
models and situations is desirable. 

Let us now analyze the evolution of the mesonic properties in connection with the possible 
restoration of symmetries. In the present work, the mesonic behavior is driven only by the 
degree of restoration of chiral symmetry in the different sectors. This does not exclude that 
other effects, not taken into account here, might influence its behavior. It can be seen in 
Fig. [2] (upper panel) that the partners (vr, a) and (r/, ao) become degenerate at almost the 
same temperature. In both models, this behavior is the signal of the effective restoration of 
chiral symmetry in the nonstrange sector. 

On the contrary, the r]' and /o masses do not show a tendency to converge in the region 
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FIG. 3: Scalar and pseudoscalar mixing angles as a function of the reduced temperature T/Tc 
for the PNJL (solid lines) and NJL (dashed lines) models; the ideal mixing angles are also shown 
(dotted lines). 



of temperatures studied, a behavior that reflects the reluctance of chiral symmetry to get 
restored in the strange sector. In fact, due to the behavior of the mixing angles, /o and rj' 
become essentially strange as the temperature increases. Moreover, as it has been shown in 
Ref. 3J], even when the dynamically broken chiral symmetry is restored in all sectors (and 
differently from what is found for nonstrange chiral partners) a sizable difference between 
the masses of these mesons survives, a fact due to the high value of the current strange quark 
mass used here (m^ = 140.7 MeV). Indeed, at high temperatures m\ ~ m?, + 4m^, thus 
explaining the observed behavior. Finally, we focus on the k meson (Fig. [2] lower panel): it 
is always an unbound, resonant state and, as the temperature increases, it tends to become 
degenerate in mass with the K meson, but at a temperature of the order of 1.5 T^ (in PNJL, 
and higher in NJL). In summary, the masses of the mesons that become less strange, a and 
r], converge, respectively, with those of the non strange, 7r° and a^. The convergence of the 
chiral partners k and K, which have a us structure, occurs at higher temperatures, and is 
probably slowed down by the small decrease of the strange quark mass, Mg. 

Concerning the axial symmetry, its effective restoration should be signaled by the van- 
ishing of the observables related to the anomaly, like the mixing angles, the gap between the 
masses of the chiral partners of the U^(l) symmetry, and the topological susceptibility. For 
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the observables so far analyzed, we notice that although in both models the anomaly effects 
exhibit a tendency to decrease, a full restoration of the axial symmetry is not achieved: the 
masses of the partners {n,!]) and (cr, oq), although getting close at high temperatures, do 
not converge and the mixing angles never reach the ideal values. This was indeed expected 
since, in the framework of the NJL model, it has been shown that only with additional 
assumptions (for example, by choosing a temperature dependent anomaly coefficient |6| or 
by using a regularization where the cutoff goes to infinity at T 7^ [34]) the restoration of 
the axial symmetry can be achieved. 

Let us now comment on the differences between the results of the two models. The new 
feature of the PNJL model is that the faster decrease of the quark condensates leads to a 
faster partial restoration of chiral symmetry. Also the analysis of the mesonic masses shows 
that a faster effective restoration of this symmetry, in the nonstrange sector, is achieved, as 
can be seen in Fig. [2l In fact, in the NJL model the effective chiral symmetry restoration 
for the nonstrange sector occurs at Tg// = L3 T^ while, in the PNJL model, Te// = 1.2 T^ 
(again we do a relative comparison between the two models: Te// and T^ are derived and 
compared for each model, respectively); for the K — n sector, the temperatures are about 1.5 
(NJL) and 2 (NJL) times the corresponding characteristic temperature. Finally, although 
the axial chiral partners do not converge, in the PNJL model, the masses of (vr, rj) become 
closer than in NJL, as well as those of (a, ao). From Fig. [3] we can also see that around 
T^ the mixing angles 6p and 63 approaches faster the ideal angle in the PNJL model than 
in the NJL one. This is an indication that, although axial symmetry is not restored in the 
range of temperatures studied, the tendency to restore this symmetry is slower in the NJL 
model. 

D. Coupling constants 

In FigHlwe plot the values of the vr, K, rj and a coupling constants. We observe a striking 
behavior at the Mott temperature for each meson: the coupling strengths approach zero for 
T -^ T^^^°'^ J32|. This is due to the fact that the polarization displays a kink singularity, 
which can also be seen in the meson masses. For the 7] and a coupling constants there is 
a second drop toward zero when the mass of these mesons approach ujg = 2Mj,. In other 
terms, these two zeros signal the entrance into the continuum for u, d quarks and s quarks. 
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7'"" /T^ 


T^iott/T? 


TZou/T? 


TLu/T? 


PNJL 


1.07 


0.95 


0.98 


1.06 


NJL 


1.08 


0.82 


0.92 


1.07 



TABLE III: Reduced Mott temperatures in the NJL and PNJL models at zero chemical potential. 

respectively. 

As already stated, the most striking difference between the NJL and PNJL models lies 
in the faster variation with the temperature of the PNJL results around any characteristic 
temperature. In particular, close to the phase transition, the NJL and PNJL calculations 
for the meson-quark coupling constants show a remarkable difference. We observe that, in 
both models, the mesons without flavor mixing (vr and K) have a higher T^j^^°'^ /T^ ratio 
than those with flavor mixing {rj and a) as it can be seen in Table IIIII It is interesting to 
note that, while T^j^^"^ /T^ for vr and K does not change appreciably from one model to 
the other, for r] and a this ratio is higher in the PNJL model, where a faster decrease of 
the mixing effects is observed. This effect, which indicates a slightly longer survival of these 
mesons as bound states, is probably driven by the faster decrease of the strange quark mass, 
observed in the PNJL model. 

The PNJL model is a quantitative step toward confinement with respect to the NJL 
quark model because the $ factor suppresses the 1- and 2- quarks Boltzmann factor at low 
temperature. The fast restoration of the Z3 symmetry ($ goes to one when temperature 
increases) producing, in a short range of temperatures, a quark thermal bath with all (1-, 
2- and 3-) quark contributions might explain the fastening of the transition. 

E. Topological susceptibility 

We found it interesting to derive the topological susceptibility x? which, in pure color 
SU(3) theory, is related to the rj' mass through the Witten-Veneziano formula 361] 



J n 



-X = M' + M' - 2Mf, 



(22) 



This observable, together with the mesonic masses and mixing angles, is strongly influ- 
enced by the anomaly: besides the degeneracy of the axial chiral partners and the recovery 
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FIG. 4: The vr, K (left panel) and 77, a (right panel) coupling constants as functions of the reduced 
temperature T/T^ for the PNJL (solid lines) and NJL (dashed lines) models. 
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FIG. 5: Topological susceptibility as a function of the reduced temperature T/Tc for the PNJL 
(left panel) and NJL (right panel) models. 



of the OZI rule (mixing angles — *> ideal values), the vanishing of x is an indication of the 
absence of the anomaly and, consequently, of the effective restoration of the Uyi(l) sym- 
metry. Lattice calculations indicate a strong decrease of the topological susceptibility with 



increasing temperature 



u 



38 



39|. 
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In Fig. [5] we show that, as expected from the previous analysis of the masses of the chiral 
partners and mixing angles, the axial symmetry is not fully restored at high temperatures, 
and the topological susceptibility is far away from being zero in both models (see Fig. \^. 
However, at T = 2T^ the PN JL topological susceptibility is reduced to about 5% of its value 
at zero temperature, while the NJL one has a slower decrease. Moreover, it is interesting to 
notice that the PNJL calculation (without any change of the parameters previously fitted, or 
extra Ansatz, like the temperature dependence of the anomaly coefficient) nicely reproduces 
the first lattice points, namely, the rather steep drop around T^, while the NJL model fails 
to do so. The faster decrease of the topological susceptibility, the tendency of the mixing 
angles to approach earlier the ideal values and the reduction of the mass gap of the axial 
chiral partners, discussed in Sec. Ill C, are consistent indications that the Polyakov loop 
leads to a faster tendency to the restoration of the axial symmetry. 

IV. SUMMARY AND CONCLUSIONS 

In this work we have explored the thermodynamical properties of the vacuum state and 
the dynamics of the scalar-pseudoscalar meson spectrum propagating in a hot medium in 
the context of a SU/(3) PNJL model. Within the framework of such a model, we have 
included the fiavor mixing and the coupling of quarks to the Polyakov loop, which in turn 
is governed by an effective potential. 

Our results indicate that the main feature of the quark masses is a faster drop around 
T^ in the PNJL model than in the NJL one. This indicates that the partial restoration of 
the chiral symmetry is more efficient and fast in the PNJL model. The mass of the strange 
quark in the PNJL model is still far from reaching the strange current quark mass even for 
high temperatures, although it exhibits a faster decrease than in the NJL model. This fact 
contributes to a faster decrease of the anomaly effects, with implications in the behavior of 
several observables. 

We have observed that, qualitatively, the behavior of mesonic masses in the PNJL model is 
similar to the corresponding one in the NJL model. However, we notice that the mixing angle 
approaches the ideal angle faster in the PNJL model. In addition, the Mott temperatures 
are different in both models, showing that the domain where mesons with flavor mixing 
are bound states is extended in the PNJL model. These results show the relevance of the 
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effects of the interplay among U^(l) anomaly, the Polyakov loop dynamics, and the partial 
restoration of the chiral symmetry at finite temperature. 

As a signal of the effective restoration of chiral symmetry in the nonstrange sector, the 
partners (tt, a) and [r], ag) become degenerate, but this occurs at a lower reduced temperature 
in the PNJL model. On the contrary, in both models the rj' and /o masses do not show 
a tendency to converge in the region of temperatures studied, an indication that chiral 
symmetry is not likely to be restored in the strange sector. 

The comparative results of the meson-quark coupling constants are also interesting. In 
particular, close to the phase transition, the NJL and PNJL calculations for the meson-quark 
coupling constants show meaningful differences. In particular the a and rj mesons exhibit a 
tendency to a slightly longer survival as bound states. 

Finally, there is a significative improvement in the PNJL model in the results concerning 
the topological susceptibility. At T = 2T^ the latter is reduced to about 5% of its value 
at zero temperature, while in the NJL model it exhibits a slower decrease. Moreover, it is 
interesting to notice that the PNJL calculation nicely reproduces the first lattice points, with 
a rather steep drop around T^, while this feature is not found in the NJL model. Although 
restoration of axial symmetry is not achieved, this behavior of the topological susceptibility 
(and of other relevant observables) , indicates that the PNJL model shows a more efficient 
mechanism for the restoration of this symmetry. 
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